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The two dimensional Heisenberg antiferromagnet on the square lattice with nearest ( Ji) and next- 
nearest (Ja) neighbor couplings is investigated in the strong frustration regime (J2/J1 > 1/2). A 
new effective field theory describing the long wavelength physics of the model is derived from the 
quantum hamiltonian. The structure of the resulting non linear sigma model allows to recover the 
known spin wave results in the collinear regime, supports the presence of an Ising phase transition 
at finite temperature and suggests the possible occurrence of a non-magnetic ground state breaking 
rotational symmetry. By means of Lanczos diagonalizations we investigate the spin system at T — 0, 
focusing our attention on the region where the collinear order parameter is strongly suppressed by 
quantum fluctuations and a transition to a non-magnetic state occurs. Correlation functions display 
a remarkable size independence and allow to identify the transition between the magnetic and non- 
magnetic region of the phase diagram. The numerical results support the presence of a non-magnetic 
phase with orientational ordering. 

PACS numbers: 75.10.Jm 



I. INTRODUCTION 

Frustrated low-dimensional spin systems are still ex- 
tensively investigated since they show rich phase dia- 
grams and exhibit unusual quantum phases. A typical 
example is the square lattice J\ — Ji model: a Heisenberg 
antiferromagnet with competing couplings (Ji, J2 > 0) 
between nearest neighbor (< >) and next-nearest neigh- 
bor (<C spins: 

U = h • + J 2 ]T Si • S k (1) 

<i.j> «i,fc> 

where are spin operators. By varying the frustra- 
tion ratio a = J2 / Ji , quantum phase transitions be- 
tween magnetically ordered and disordered phases can 
take place at T — 0. The interest in the 2D frustrated 
antifcrromagnets has deep theoretical motivations in the 
characterization of disordered spin liquids (or Bose liq- 
uids). Moreover, several magnetic materials have been 
synthesized nowadays where frustration plays a dominant 
rolei. In particular, the specific interest on the frustrated 
Heisenberg model on a square lattice raised with the dis- 
covery of vanadate compounds, whose magnetic behavior 
is likely to be described by the J\ — J2 hamiltonian^. 

In the classical limit (S — ► 00) at weak frustration, 
the ground state (GS) has conventional Neel order with 
magnetic wave vector Q = (tt, n) for a < 0.5. Above this 
threshold, the two sublattices are antiferromagnetically 
ordered but remain free to rotate with respect to each 
other and the GS manifold has an 0(3) x 0(3) degener- 
acy, larger than expected on the basis of the 0(3) symme- 
try of the hamiltonian. Weak quantum fluctuations, via 
the order by disorder mechanism, can be included by use 
of spin wave theory, and are shown to select a collinear 
ordered state with magnetic wave vector Q = (jr, 0) or 
(0, 7r), reducing the GS degeneracy to 0(3) xZ 2 . In fact, 
the collinear state breaks both the 0(3) spin rotational 
invariance of the Heisenberg hamiltonian and the 7r/2 



rotational symmetry of the square lattice {Z2). In the 
S — > 00 limit, according to spin wave theory^, the criti- 
cal coupling a c = 1/2 marks the first order transition be- 
tween the collinear and the Neel phase. However, when 
quantum fluctuations are taken into account beyond per- 
turbation theory, (at least) one intermediate phase is ex- 
pected to separate the two magnetically ordered phasesi. 

Nevertheless, on the basis of a very recent pertur- 
bative numerical renormalization group analysis no ev- 
idence for an intermediate phase has been found and it 
has been proposed that a direct and unexpected second 
order phase transition may occur at the classical critical 
poinl4. 

The aim of this work is to investigate the quantum phase 
transition occurring at T = as the frustration ratio a 
is decreased, i.e. when the collinear order is suppressed 
by quantum fluctuations. To this purpose, the effective 
long-wavelength action of the two dimensional quantum 
model in the regime of strong frustration is mapped into a 
non-linear sigma model (NLSM) in 2+1 dimension. This 
is the first time that the mapping, originally proposed by 
Haldane for a Heisenberg chain^, is properly generalized 
to the strongly frustrated J\ — J2 model. On the basis 
of the symmetries of the NLSM action, an Ising phase, 
breaking the J rotational symmetry of the square lat- 
tice and preserving the 0(3) spin rotational invariance, 
may exist both at low and zero temperature. While the 
stability of this phase at finite temperature has been pre- 
viously investigatedSiS^, it has never been explicitly ex- 
amined if such an Ising state is stable at zero tempera- 
ture for some values of the frustration ratio lower than 
the one corresponding to the onset of the collinear or- 
der. Lanczos diagonalizations support this scenario via a 
careful examination of both the excitation spectrum and 
the correlation functions: A non-magnetic valence bond 
nematic phase with orientational ordering is the most fa- 
vorable candidate as a ground state in a portion of the 
phase diagram. 
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II. NON LINEAR SIGMA MODEL FOR THE 
COLLINEAR PHASE 

Among the various theoretical approaches adopted for 
the J\ — J 2 model in the regime of weak frustration, the 
NLSM method is particularly suitable for the study of 
the phase transition between a magnetically ordered state 
and a disordered phase. The 2D frustrated Heisenberg 
antifcrromagnet with a < a c is mapped to a 0(3) NLSM 
in D = 2 + 1 dimensioni2*ii, which indeed shows a second 
order (quantum) phase transition to a non-magnetic state 
at T = 0. When frustration is strong and the GS is 
collincar, the mapping to a 3D classical model is still 
possible, but the effective long- wavelength action is no 
longer a conventional 0(3) NLSM. 

Here we generalize the original mapping proposed 
by Haldane for the microscopic derivation of the 
long- wavelength, low-energy effective theory of one- 
dimensional quantum antiferromagnets in the Neel 
phased, to the J\ — J 2 model in the strong frustration 
regime, where collincar order is expected. By using the 
Trotter formula and a coherent state basis in the spin 
Hilbert space, the partition function of the system is writ- 
ten in a path integral representation as: 



T>Cl exp 



[sp]) 



(2) 



where Oj(t) is a classical 0(3) vector field defined on each 
lattice site, normalized to Clf(t) = 1 and S is the action 
at a temperature 

rP 

S[Cl(T)} = -iSy2 UJ Pi( T )]+ drH[fi(r)] (3) 

Jo 

The first (purely imaginary) term is the Berry phased and 
H[f2] is the expectation value of the hamiltonian operator 
on the coherent states basis: 



In the regime of strong frustration it is convenient to 
separate the lattice in the two sublattices which will be 
labeled + and — respectively Following Haldane^, we 
split the spatially oscillating spin state fli on the sub- 
lattice + (— ) as the sum of two orthogonal smooth vec- 
tor fields, describing the local Neel order n + (n_) and 
the transverse fluctuations L + (L_), satisfying the con- 
straints = 1 (n_ = 1) and n + • L + = (n_ • L = 0). 
In order to carry out the splitting between the uniform 
and staggered fluctuations keeping the right number of 
independent variables, we proceed by partitioning the 
lattice in plaquettes as shown in Fig. ^ Then we define 
the Neel fields (n+, n_) and the associated fluctuations 
(L_|_, L_) in the center of each plaquette as follows: 



Mr) 2 , Mr) 



S 



s 



(5) 




FIG. 1: Lattice of the Ji — J2 model. The lattice spacing 
is a. Grey and black circles denote + and — sublattice sites 
respectively, while white circles the center of each plaquette, 
where the fields n+, n_ and L+, L_ are defined. 
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where 77 = ±1, a is the spacing of the original lattice and 
r = (x, y) is the coordinate of the plaquette centers which 
define a square superlatticc of spacing 2a. In the contin- 
uum limit, to second order in space and time derivatives 
and keeping the lowest order in 1/5, the partition func- 
tion is written as 



Z = J Vh± X>L± 5 (n+ • L+)<5 (n_ • L_ 



(7) 



where 
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d 2 x I drC 





and the lagrangian is 

C = L T AL - B T L + Ki 



K 2 



(8) 



T>h± and Vh± mean integration over the Neel fields and 
the fluctuations in both the sublattices and the six com- 
ponent array L T = (L + ,L_) gathers the (three) com- 
ponents of the two fluctuation fields L" and L". The 
scalars K are defined as 

Ki = -S 2 Ji (d y h + ■ d y n_ - d x n + ■ d x n_) 

K 2 = S 2 h [( VA ^) 2 + k ■ d y n k ] (9) 

k=± 

while the 6x6 matrix A is written in block form as 



. 1_ ( 2J 2 h 

a 2 \ Ji 2J 2 

and the array B T = (B + ,B_) is given by: 

B+ = — (n+ x d T n+) + 7 _ + 74 

zgt a a 

B_ = —2 (n_ x S T n_) + 7+ + 7 _ 

2a^ a a 



(10) 
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with 



7+ 



7- 



&.n_ 



To this order in derivatives the Berry phase contribu- 
tions iw[n+(r)] and iw[n_(r)] identically vanish^. Per- 
forming the gaussian integration of the partition function 
with respect to fluctuations needs some care due to the 
constraints n + • L + = and n • L = which limit the 
integration to the transverse components of the spin fluc- 
tuations and L^. In order to release the constraint, 
we first multiply the partition function by a constant fac- 
tor F (which does not affect the physical properties of the 
model) written as a gaussian integral over two auxiliary 
scalar fields v+ and w_: 



F = J Vv±e S d2x Xf 



dr (L 11 ) AL" 



where L|! = v + h + (L_ = w_n_). By defining the vector 



L" (L. 



Li) , the integral over 



fluctuations X in Eq. J7J) can be written as: 
J = [ VL ± e~f d2x J» ^[( lX ) t al^+( L iI)^ L h -b*V 



where L;J: is explicitly given by LzJ: = L + — (L + • n + ) n + 
( = L — (L_ • n_)n_). As a result, the partition 
function is written as an unconstrained integral over the 
fields L± and h±. The lagrangian maintains the same 
formal structure JHJ with a modified matrix A and array 
B: 



A 



a/3 



A 



/ k 



rap 



X - n%nl + 2(n t ■ n k )nfn% 



B; 



Bf-(B i -n i )n? 



(11) 



where Latin indices i, k = ± identify the sublattice and 
Greek superscripts run over the three spin components. 
In practice, the gaussian integration has been performed 
by diagonalizing the matrix A on its eigenvector basis 
{uf (I)} i with I = 1, . . . , 6 so that the effective lagrangian 
density is expressed in terms of the eigenvalues {A/} ; 

C = K 1 +K 2 -lj2 X i 1 bf 6/ = ^Bf<(0 (12) 



After tedious but straightforward calculations, the result- 
ing effective lagrangian is written as 



£[n + ,n_] = £ s +£j 



(13) 



with 
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i=± 
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2 
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2a 2 D 



(<9 y n+ • d y h- - d x n+ ■ d x n-) 
(n+ • 9 T n_) (d T n + ■ n_) + 



Ji 2 J 2 r 

8a 2 D 
Ji 
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Where D* = 4 J 2 2 - a 2 h 2 , D = D*{AJ 2 2 - J t 2 ) and 
(7 = n + ■ fi_ is the "Ising order parameter" . The la- 
grangian density C preserves the 0(3) symmetry of the 
microscopic hamiltonian Q but breaks the invariance of 
the model under ir/2 lattice rotation due to the adopted 
coarse graining procedure^. However, the original rota- 
tional invariance reflects in the additional Z 2 symmetry 
of our lagrangian under the simultaneous action of tt/2 
rotation and inversion of one field: (x, y) — > (— y, x) and 
n+ — > — n_|_. The global symmetry group of the resulting 
field theory is therefore 0(3) x Z 2 - 

In order to check the correctness of the effective la- 
grangian density C, we compare the dispersion relation 
arising from a saddle point evaluation of the partition 
function to the known spin wave results. The saddle 
point configuration corresponds to the minimum of the 
action, i.e. to a homogeneous and static configuration 
characterized by two independent unit vectors and 
n°_ which describe one of the degenerate classical ground 
states. By expanding up to second order in fluctuations 
n + = n , + 5ii + (n_ = n°_ + <$n_) and performing the 
gaussian integration we get precisely the same result of 
spin wave theory in the long wavelength limit^ii4: The 
low energy excitations are described by four branches la- 
beled by A = ±1 and rj = ±1 whose dispersions are given 
by 



0J 2 {k) = 4S 2 {k 2 x {2J 2 +r 1 J 1 )(2J 2 + ? ] J 1 cos0) + 

" 2 <^ T - 1 AJi) (2J 2 - AJiCOS0)} (14) 



k z y (2 J 2 



where 9 is the angle between the two Neel fields: cos 8 = 
n 1 -n_ . Accordingly, perturbation theory on the effective 
action shows that the lowest free energy is attained when 
the staggered magnetizations of the two sublattices n*]_ 
and n_ are either parallel or antiparallel: The collinear 
order is stabilized by quantum fluctuations. This picture 
is believed to be correct in the limit of large frustration 
ratio a. By lowering a, fluctuations are enhanced and 
a quantum phase transition is expected to occur at zero 
temperature before reaching the (classical) limiting value 
a = 0.5. On the basis of the symmetries of the effective 
model, an intermediate regime characterized by vanishing 
staggered magnetization < n + > = < n_ > = and pos- 
sibly a finite Ising order parameter < a >= n + • n_ 7^ 
may be present. This hypothetical state would break the 
rotational symmetry of the lattice preserving the SU{2) 
spin symmetry. Due to the coarse graining carried out 
in the microscopic hamiltonian, such a state may cor- 
respond either to a translationally invariant " valence 
bond nematic" phase where valence bonds display ori- 
entational ordering or to a valence bond crystal (VBC) 
which breaks the translational symmetry. 

A low energy, long wavelength effective field theory for 
collinear antiferromagnets was put forward in a seminal 
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paper by Chandra. Coleman and Larking. At variance 
with our approach, in that work quantum fluctuations 
were integrated out by use of a perturbative spin-wave 
approximation and the analysis was limited to the finite 
temperature domain. In this way, CCL predicted a finite 
temperature Ising transition in two dimensions: Even if 
at T = both sublatticcs display finite staggered mag- 
netization in zero field, long-range Neel order disappears 
at any finite temperature, since thermal fluctuations are 
known to restore the continuous symmetries in two spa- 
tial dimensions^. CCL suggested that the previously 
defined Ising order parameter a may preserve long-range 
order up to a non-zero critical temperature. Recently the 
presence of such a finite temperature transition has been 
the subject of several investigation:,-------. 

The CCL approach, however, cannot be directly ap- 
plied to the study of the zero temperature limit be- 
cause the effects of quantum fluctuations are consid- 
ered only within perturbation theory. In order to clarify 
the relationship between our effective field theory and 
the CCL approach, it is convenient to specialize the 
lagrangian £[n + ,n_] to a class of field configurations 
of the form n(r, t)+ = n(r)9_ + <$n(r, f)+ (n(r, t)_ = 
n(r)1 + <5n(r,i)_) characterized by weak fluctuations 
(<5n(r, i)) on top of a slowly varying time independent 
(i.e. classical) configuration (n(r)). Expanding C to sec- 
ond order in the fluctuations and performing the gaus- 
sian integration, by use of the result lfl"4"|l , we recover the 
CCL result. This approach can be justified at large frus- 
tration ratios a (i.e. deep in the collinear regime) where 
quantum fluctuations are not able to severely affect the 
classical ground state. However, in order to analyze the 
quantum transition between collinear order and a disor- 
dered phase at T = 0, it is necessary to take into account 
the effects of quantum fluctuations beyond the spin wave 
approximation, i.e. we have to study the full effective 
lagrangian (II 31) . 



III. LANCZOS DIAGONALIZATIONS 

By means of Lanczos diagonalizations (LD) we try to 
clarify if an Ising state can be stabilized when the long 
range collinear order fades away, as suggested by the field 
theoretical approach. 

Lanczos diagonalizations have been performed on the 
4x4 and 6x6 square clusters for spin S = 1/2. By 
means of LD we obtain the energy spectrum, providing 
indications of possible changes in the nature of the GS 
which occur by increasing the frustration ratioi&ii. In 
the collinear phase, due to the spatial symmetry break- 
ing, four classes (towers) of states with different spatial 
symmetries are expected to become degenerate: the low- 
est representative of these classes are an s-wave and a 
d-wave singlet at momentum (0, 0) and two triplets at 
momenta (0, w) and (ix, 0). The low energy spectrum 
as a function of the frustration ratio is shown in FigEJ 
While the behavior of the model for a ~ 0.5 has been the 
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FIG. 2: Low energy states referenced to the GS as a function 
of the frustration ratio a = J2/J1 in the 4x4 (upper panel) 
and 6x6 (lower panel) cluster. 

subject of several investigations^ and is still a debated 
problem, here we will concentrate on the GS properties 
in the regime of larger frustration a > 0.6. A clear ten- 
dency to break the rotational symmetry is suggested by 
the quasi degeneracy of the s-wave and d-wave singlets 
which actually cross each other in the 6x6 cluster. How- 
ever, this does not rule out the possible occurrence of the 
collinear phase down to a — 0.6, since the energy gap 
of the triplet at (0, it) is shown to decrease with the lat- 
tice size. The (ir, 0) singlet gap also decreases with size 
although it remains considerably larger than the lowest 
singlet gap and in fact comparable to the (n, ir) triplet 
gap which is believed to be finite in the thermodynamic 
limit for a > 0.5. The (tt, 7t) singlet (not shown in figure) 
is much higher in energy Therefore, from the ordering 
of the low energy states we may conclude that: i) rota- 
tional symmetry is broken for a > 0.6; ii) triplet states 
are gapped for a < 0.62; Hi) the columnar VBC phase 
is unlikely to occur, at least for a > 0.62 and iv) other 
VBC phases, like a plaquette state-----------, are not com- 
patible with the observed ordering of levels in the energy 
spectrum. 

The quasi degeneracy between the s-wave (\s >) and d- 
wave (\d >) singlets in the extended range 0.6 < a < 0.7 
allows to consider the two real linear combinations of 
these states as good representations of the symmetry bro- 
ken phases which are physically realized in the thermo- 
dynamic limit. This is particularly convenient since in 
a non rotationally invariant state, like (|s > ±\d >)/v2, 
the rotational order parameter Or = Sr • Sr+g — Sr • Sr+x , 
(where x and y are the two primitive vectors of the lat- 
tice) may acquire a non-zero value. If the two singlets are 
degenerate and < O r > remains finite in the thermody- 
namic limit, rotational symmetry breaking occurs in the 
model2--. In the thermodynamic limit, this procedure 
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FIG. 3: Properties of the symmetry broken state: Ising order 
parameter O (triangles); Probability to find next neighbors 
singlets in a given direction; P M =< 1 — §(Sr + Sr+p) 2 >• 
Empty dots: P x , Full dots: P y . Dashed line: P x for indepen- 
dent Heisenberg chains. 




berg chain. In essence, the whole system seems to behave 
as a collection of spin chains weakly coupled in the trans- 
verse direction. As a grows, P x keeps almost constant, 
while P y decreases, leaving room to the formation of ver- 
tical triplets: The system is moving toward the collinear 
phase. 

In order to better characterize the phase diagram we 
also investigated the behavior of spin correlations. The 
Fourier transform S(k) of < SqS* > along a closed path 
in the Brillouin zone has been calculated both for the 
s-wave and d-wave singlet and the results are shown in 
Fig. The close similarity between the spin correlations 
in the two states indeed confirms that they both contain 
the same physics. At J2 = 0.55 J\ S(k) exhibits a peak 
at momentum (7r,7r) suggesting that the dominant (short 
range) correlations are still antiferromagnetio^. In the 
range 0.60 < a < 0.62 S(k) is remarkably flat and does 
not show significant size dependence. Instead it seems 
that for larger frustration a the system is going to sus- 
tain a transition to a collinear phase. Finally, at a = 0.7 
the (71", 0) peak in fact grows quite substantially with the 
size, signaling the onset of the collinear order parameter. 
A remarkable common feature of the LD results is the 
collapse of the 4x4 and 6x6 data on the same smooth 
curve except (possibly) at a single wave- vector, which 
identifies the dominant periodicity in the spin correla- 
tions. The LD data then allow for an accurate evaluation 
of the full momentum dependence of the magnetic struc- 
ture factor: we performed a fit of the numerical data (also 
shown in Fig. 0J) with a parameterized form inspired by 
the spin wave theory results^. The chosen function rep- 
resents quite accurately the numerical data except at the 
single wave- vector where the order parameter sets in and 
a singular contribution develops in the thermodynamic 
limit. 



FIG. 4: Magnetic structure factor along a closed path in the 
Brillouin zone from T = (0,0) to M = (0,tt), to X = (n, vr) 
and back to T. Full (empty) triangles: S(k) evaluated on the 
lowest s-wave (d-wave) singlet for a 4 x 4 cluster; Full (empty) 
dots: S(k) evaluated on the lowest s-wave (d-wave) singlet for 
a 6 x 6 cluster 

would be fully equivalent to the usual way to evaluate 
order parameters in terms of the asymptotic behavior of 
correlation functions < O^Oq >. However, in small clus- 
ters, we believe that our approach is less affected by finite 
size effects. The two (quasi) degenerate states have van- 
ishing momentum and then the order parameter < Or > 
is translationally invariant i.e. independent of r. The 
numerical results arc displayed in Fig. [3] together with 
the probability to find next neighbor singlets both on 
horizontal P x and vertical P y bonds. 

A strong anisotropy is present: The order parameter 
is large and the singlet probability strongly differs in the 
two spatial directions. Remarkably, around a ~ 0.60, 
P y ~ 0.25 is compatible with a disordered configuration, 
while P x gets close to the limit characterizing the Heisen- 



IV. CONCLUSIONS 

In this paper we derived, for the first time, the long 
wavelength, low energy effective field theory describing 
quantum and thermal fluctuations in the collinear phase 
of frustrated 2D antifcrromagncts. The resulting NLSM 
is written in terms of two fields describing the local Neel 
order parameter of the two sublattices and is invariant 
under the 0(3) x Z2 symmetry group. On the basis of 
this formalism we are led to predict the possible occur- 
rence of a non-magnetic ground state breaking rotational 
symmetry for suitable values of the frustration ratio. In 
order to investigate this possibility, we also performed 
Lanczos diagonalizations. By a careful inspection of the 
numerical results we found evidence for the occurrence 
of the predicted valence bond nematic ground state in 
a region around a ~ 0.6. The evaluation of the mag- 
netic structure factor in small clusters also showed that 
the short range spin correlations of the J\ — J2 model 
are remarkably size independent. This observation may 
be very useful in the interpretation of accurate neutron 
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scattering data on frustrated 2D antifcrromagncts. 

According to spin wave theory the J\ — J2 model dis- 
plays a first order phase transition at a c = 0.5 between a 
Neel and a collinearly ordered region in the classical limit 
(S — > 00). When quantum fluctuations are taken into 
account an intermediate SU (2) invariant phase is stabi- 
lized. Many different candidates have been proposed as 
possible ground states in this region: gapped or gapless 
spin liquids 22 ! 24 , VBC's with columnar 16 i 19 i 20 i 25 or pla- 
quette pattern s 11 ! 18 ! 19 . 

At any finite temperature, the continuous spin rota- 
tional symmetry, broken in the collinear phase, is re- 
stored and only the breaking of the symmetry corre- 
sponding to the order parameter a = n + ■ n_ can in prin- 
ciple survive up to a finite critical temperature defining 
a phase transition which lies in the 2D fsing universality 
clasA The hypothesis that as T — > the transition line 
ends in a point different from the one corresponding to 
the onset of the collinear order can not be excluded a pri- 
ori. In such a case, the ground state in a portion of the 
intermediate 577(2) invariant region may be a valence 
bond nematic phase with some orientational ordering. 
This possible scenario is consistent with the suggestions 
of the analytical results based on the NLSM action and 
has been confirmed by a LD analysis. 



Because of the quasi degeneracy in the energy spec- 
trum in a region around a ~ 0.60, a state characterized 
by an orientational symmetry breaking is very likely to 
occur, ruling out the fully symmetric spin liquid. Simi- 
larly, the ordering of excited states is not compatible with 
the plaquette VBC. A careful analysis of the low energy 
spectrum and of the correlation functions suggests that 
a zero temperature transition takes place at a c ~ 0.62. 
The transition separates the large a collinear phase and 
an intermediate gapped regime breaking the it/2 rota- 
tional symmetry of the lattice. Thus on the basis of our 
investigations we argue that such a phase may be conve- 
niently thought of as "nematic" ordering of valence bonds 
and anticipates an isotropic spin liquid (or a VBC) which 
is likely to occur at lower a. 

A direct transition between the collinear phase and a 
VBC should be of the first order, the two phases having 
different order parameters. Instead the transition to the 
Ising phase may be of second order being related just to 
the vanishing of the sublattice staggered magnetization. 
Our analysis is not able to discriminate between second 
order and weakly first order transition: Monte Carlo sim- 
ulations of the NLSM action derived here will be helpful 
to supplement LD data in order to clarify this issue. 



1 See for instance: Misguich and C. Lhuillier in "Frustrated 
spin systems", H. T. Diep editor, World-Scientific (2005) 
and cond-mat/0310405 and Refs. therein. 

2 R. Melzi et ai, Phys. Rev. Lett. 85, 1318 (2000); Phys. 
Rev. B 64, 024409 (2001). 

3 P. Chandra a nd B. Doucot, P hys. Rev. B 38, R9335 (1988). 

4 S. Moukouri, cond-mat/0504306 (2005). 

5 F. D. M. Haldane, Phys. Rev. Lett. 50, 1153 (1983). 

6 P. Chandra, P. Coleman, A. I. Larkin, Phys. Rev. Lett. 
64, 88 (1990). 

7 R.R.P. Singh, W. Zheng, J. Oitmaa, O.P Sushkov, C.J. 
Hamer, Phys. Rev. Lett. 91, 017201 (2003). 

8 C. Weber et ai, Phys. Rev. Lett. 91, 177202 (2003). 

9 L. Capriotti , A. Fubini, T. Roscilde, V. Tognetti, Phys. 
Rev. Lett. 92, 157202 (2004). 

S. Chakravarty, B. I. Halperin and D. R. Nelson, Phys. 
Rev. B 39 2344 (1989). 

11 K. Takano, Y. Kito, Y. Ono and K. Sano, Phys. Rev. Lett. 
91, 197202 (2003). 

12 This is true within the 0(3) NLSM approach where it is 
assumed that the Neel field is well defined at all points in 
(2+1) dimensional space-time. However if this assumption 
is relaxed to allow "hedgehog" (spatial) singularities of the 
Neel field then a dependence on the quantized value of the 
spin S appears and the Berry term may contribute to the 
Boltzmann weight of the configuration—. 



F. D. M. Haldane, Phys. Rev. Lett. 61, 1029 (1988). 

A. Moreo, E. Dagotto, T. Jolicoeur and J. Riera, Phys. 

Rev. B 42, 6283 (1990). 

N. D. Mermin, H. Wagner Phys. Rev. Lett. 17, 1133 
(1966). 

D. Poilblanc, E. Gagliano, S. Bacci and E. Dagotto, Phys. 
Rev. B 43, 10970 (1991). 

T. Einarsson and H. J. Schulz, Phys. Rev. B 51, R6151 
(1994); H. J. Schulz and T. A. L. Ziman, J. Phys. I 6, 675 
(1996). 

M. E. Zhitomirsky, K. Ueda, Phys. Rev. B 54, 9007 (1996). 
O. P. Sushkov, J. Oitmaa and Z. Weihong, Phys. Rev. B 
63, 104420 (2001). 

R. R. P. Singh, Z. Weihong, C. J. Hamer, J. Oitmaa, Phys. 
Rev. B 60, 7278 (1999). 

These two conditions together imply that the susceptibility 
associated to the order parameter O diverges. 
L. Capriotti, F. Becca, A. Parola and S. Sorella, Phys. Rev. 
Lett. 87, 097201 (2001). 

The fitting function is the square root of the ratio between 
two first order polynomials in 7(k) = cos k x + cos k y and 
7'(k) = cos k x cos k v . 

F. Figueirido et al. Phys. Rev. B 41, 4619 (1990). 

V. N. Kotov, J. Oitmaa, O. P. Sushkov and Z. Weihong, 

Phys. Rev. B 60, 14613 (1999). 



